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Abstract 

> 

We present three-dimensional, non-axisymmetric distorted black hole initial 
data which generalizes the axisymmetric, distorted, non-rotating |jj and ro- 
tating (2| single black hole data developed by Bernstein, Brandt, and Seidel. 
^ ■ These initial data should be useful for studying the dynamics of fully 3D, 

distorted black holes, such as those created by the spiraling coalescence of 

O ! two black holes. We describe the mathematical construction of several fam- 

ilies of such data sets, and show how to construct numerical solutions. We 

^Jy survey quantities associated with the numerically constructed solutions, such 

as ADM masses, apparent horizons, measurements of the horizon distortion, 

•<-h . and the maximum possible radiation loss (MRL). 

■ 
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I. INTRODUCTION 



Black hole studies have received significant attention in numerical relativity over the past 
several years, as the computers needed to solve the Einstein equations grow in power. The 
need for computer generated templates of gravitational waves has become more pressing 
since gravitational wave detectors, such as LIGO and GEO600, are nearing completion and 
should begin taking data in about a year Binary black hole coalescence events are 

considered to be prime candidates for the first detection of gravitational waves ||. 

Several theoretical approaches have been developed for treating these systems. So far, 
the post-Newtonian approximation (PN) has provided a good understanding of the early 
slow adiabatic inspiral, or "far-limit" , phase of these systems [|6|,[7|,|8|,|9|JT0|| . Similarly, for the 
final moments, when the black holes are close enough to each other to sit inside a common 
gravitational well, one can successfully apply the "close limit" approximation (CL) ||11|| , 
which effectively describes the whole system as a perturbation of a single black hole which 
rapidly "rings-down" to stationarity. Before this last stage, though, when the black holes 
are still close to the innermost stable circular orbit, the orbital dynamics are expected to 
yield to a plunge and coalescence. No approximation method can be applied in this highly 
nonlinear phase and it is generally expected that one can only treat the system by a full 
numerical integration of Einstein's equations. 

Ideally, one would like to start a full nonlinear integration of the Einstein equations 
with initial data that correspond to black holes in the early adiabatic inspiral phase. Such 
data would be straight-forward to compute [O], and would correspond well to a realizable 



astrophysical situation. Unfortunately, there are several technical difficulties involved in nu- 
merically evolving such a system all the way through coalescence to ringdown. Accordingly, 
in order to develop and test fully nonlinear numerical codes while the difficulties are being 
resolved, the strategy has been to evolve initial data that can be used as a model of the 
merger event. For example, initial data corresponding to holes which are initially very close 
together have been evolved p3|JT^JT5[] , as have data sets which correspond to a single very 
distorted hole 
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In this paper we describe initial data sets which correspond to three-dimensional distorted 
black holes, both with rotation and without. The formalism employed extends the axisym- 
metric work done by Bernstein, Brandt, and Seidel Although these data sets and their 
evolution have already been described briefly in several places p^8| , p!9| , |20| , |21| , p2| , [23| , p^| , in this 
paper, we provide a more complete description of them. 

The paper is organized as follows. In Section |I| we describe the mathematical setup and 
numerical construction of the data sets. In Section [III] we give results of parameter studies 
of their properties, such as ADM masses, apparent horizon positions, distortion measures, 
and maximal radiation loss. Then, we conclude our results in Section ITV. 



II. CONSTRUCTING DISTORTED 3D BLACK HOLES 

The standard approach to constructing initial data for Einstein's equations is to consider 
a 3+1 split (3 space and 1 time) of the full four- dimensional theory. In this way the equations 
divide naturally into two classes: four constraint equations, and six evolution equations. The 
procedure for performing the split is described extensively in many places, notably in the 
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review article by York |[25|| . The constraints, which are elliptic in nature, relate the 3-metric 
jij and extrinsic curvature at any coordinate time t. In our case, we solve the constraints 
to obtain data sets for some initial time labeled by t = 0. The constraints can be further 
split into the Hamiltonian constraint and the momentum constraints. Here, we consider 
vacuum spacetimes, for which the Hamiltonian constraint is 

R + (trK) 2 - K ij Kij = 0, (1) 

and the three momentum constraints are 

Di(K ij - fHvK) = 0, (2) 

where R and Di are the scalar curvature and covariant derivative associated with the 3- 
metric 7y. 

To solve the constraints, it is common to use York's conformal decomposition 



method |[25|| . This method starts by factoring out a function ip, known as the conformal 



factor, from the 3-metric and extrinsic curvature tensor components in the following way: 

7*j = (3) 
K i} = ij~ 2 K ir (4) 

If we use this decomposition, and restrict ourselves to initial slices with vanishing trK = 
j^Kij, the constraint equations take the following form: 

A^ = ^R-^- 7 K i ^k ij (5) 

bik ij = dik i] + f i k k kj + f\ k k ik = o, (6) 

where A is the Laplacian, R the scalar curvature, and T\ k the Christoffel symbol associated 
with the conformal 3-metric 7^-. Note that the conformal factor does not appear in the 
momentum constraint equations. This allows us to use the following procedure to derive 
initial data. First, one specifies the conformal 3-metric %j freely. Then, the momentum 
constraint equations are solved for the conformal extrinsic curvature. Finally, one solves the 
Hamiltonian constraint for the conformal factor. This is the procedure we use to create the 
3D black hole initial data sets described below. 



A. Distorted non-rotating black hole 

As part of his thesis work, Bernstein studied initial data sets corresponding to single 
black holes which were non-rotating, distorted, and axisymmetric Q. For these non- rotating 
black holes, the extrinsic curvature was taken to vanish. In this case, the vacuum momentum 
constraint equations are satisfied identically, and the Hamiltonian constraint reduces to 

= ~v>A (7) 

The form of the conformal 3-metric was that used by Brill in his study of pure gravitational 



wave spacetimes | 26|| . Using spherical-polar coordinates, one can write the 3-metric as 
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dl 2 = V 4 [e 2q {dr 2 + r 2 d0 2 ) + r 2 sin 2 Odcj) 2 } , (8) 

where q is the Brill "packet" which takes some functional form. Using this ansatz with (|7|) 
leads to an elliptic equation for if> which must be solved numerically. The inner boundary of 
our computational domain is the throat of the black hole, and there we apply an isometry 
condition on ip which maps the region exterior to the throat to another asymptotically flat 
region interior to it. At the outer boundary, we apply the Robin condition which requires 
ip ~ 0(r _1 ) as r — ► oo is used as the outer boundary condition. The choice of q = produces 
the spherically symmetric Schwarzschild solution in isotropic coordinates. In Bernstein's 
work, a logarithmic radial coordinate 77 is used, related to the isotropic radial coordinate r 
by 

where M is the mass of the Schwarzschild black hole that results from setting q — 0. In 
this coordinate system, the 3-metric is 

di 2 = ^VW + de 2 ) + sin 2 edf], (io) 

and the Schwarzschild solution is 



# = V^cosh(|). (11) 

Note that the conformal factor in the rj— coordinate system, if), differs from that in the 
isotropic radial coordinate system by a factor of r 1 / 2 . 

The Hamiltonian constraint in this coordinate system is 

d 2 if) d 2 ip n dip 1 ~,d 2 q d 2 q lN ... 

W + W + cot °ae = + df 2 ~ iy (12) 

The throat of the black hole is located at r = Mo/2, or 77 = 0. The isometry condition 
across the throat in the 77 coordinate system takes the particular simple form 7^(77) = 
±7y(— 77), which for the conformal factor becomes ip(f]) = ip{~v)- At the outer boundary 
we require that if) have the same behavior as the spherically symmetric solution. 

The Brill wave function q was chosen to have the following form, which obeys the bound- 
ary conditions but is otherwise arbitrary. Specifically Brill showed that in order for the mass 
of the hypersurface to be well defined, q must vanish on the axis and decrease radially at 
least as fast as e~ 2r] . In this paper, we choose q to be of the form used by Bernstein: 



q(r),9)=Q sm n 9 



' v+vo ^ 2 ^ y-vg ^ 2 



(13) 



Roughly speaking, Q is the amplitude of Brill wave, rjo its radial location, and a its width. 
Regularity along the axis requires that the exponent n must be even. We characterize initial 
data sets by the parameters (Qq, rj , a, n). 

Note that an even exponent n results in data that have equatorial plane symmetry, so 
that only 6 in the range < 9 < n/2 (or equivalently, tt/2 < 6 < 71) need to be considered. 
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We chose to extend Bernstein's work to 3D by multiplying the Brill wave function q by a 
factor which has azimuthal dependence. The particular form of this factor was 1 + c cos 2 4>, 
giving the g-function 



q(r], 9,<f>) = Q sm n e 



i v+no \2 f n-vp \2 

G I G 



ccos 



(14) 



Here, one can see that the above axisymmetric case is recovered by setting c = 0. Now the 
Hamiltonian constraint takes on the more complicated form 

9 2 1p C^V cgc 2 ee 2q&± , 2 cgc 2 Q e 2q^l^_ = 

drf d6 2 86 dcj) 2 d(f> d(f> 

2 



d2( l . . o 2q 2n d2 1 , o 2q 2 n ( 9 Q V 

— + — + 2e^csc ^ + 3e«csc ^- j 



(15) 



In this case, we characterize our initial data sets by the parameters (Q , r/ , a, c, n). 

Notice that we still have equatorial plane symmetry with this q. We also have sufficient 
symmetry in the azimuthal angle <fi to restrict ourselves to the range < <fi < tt/2. This 
allows us to perform computations in only one octant, resulting in a great savings in compu- 
tational time. However, for consistency with the non-axisymmetric rotating cases in which 
one cannot apply any symmetries, we solve it on a full grid domain. Note also that although 
setting n = 2 is valid in axisymmetry, when c ^ 0, n must be at least 4 for the right hand 
side of fllBD to be regular. Using the same boundary conditions as above, we can solve fllBl) 
for the conformal factor ip numerically on this spherical grid. 

In order to evolve these initial data, we interpolate the conformal factor ip and its deriva- 
tives onto a Cartesian grid. As a test of our initial data solver, we compute the residual of 
the Hamiltonian constraint, H, on this Cartesian grid at various resolutions. We expect to 
see second order convergence in the grid spacing. 

Fig. ^ shows the logarithm of the residual of the Hamiltonian constraint, H, on the x-axis 
at three different resolutions for the data set (Q , rjo, (J, c, n) = (—0.5, 0, 1, 1, 4). The values of 
H are rescaled so that they would coincide for perfect second order convergence. We obtain 
near second order convergence, although some error is seen for the coarsest resolution. 

This family of initial data sets of isometric embedding |27| , and their evolution as pertur- 



bations of the Schwarzschild black hole [|20|j2TJ| have been studied. In Section |TJ, we study 
various physical properties of these black holes as a function of the parameters of this initial 
data. 



B. Distorted rotating black hole 



Brandt and Seidel extended the above non-rotating axisymmetric initial data sets of 
Bernstein to include rotation P JT^j28|] . In the same way that Bernstein's data sets correspond 
to a Schwarzschild black hole surrounded by a gravitational wave, their data sets correspond 
to either a Kerr or Bowen and York black hole surrounded by a gravitational wave. For the 
rotating case, the 3-metric takes the form 



dV 



?2{q-q ) 



(drf + dO 2 ) + sinfl 2 



(16) 



5 








_16 I ■ 1 ■ 1 ■ 1 ■ 1 

-4 -2 2 4 

x/M 

FIG. 1. We show the convergence of the residual of the Hamiltonian constraint for distorted 
non-rotating black hole initial data with parameter set (Qq, rjo, a, c, n) = (—0.5,0,1,1,4). H is 
plotted for Cartesian grid sizes (resolutions) 259 3 (0.024M), 131 3 (0.047M), and 67 3 (0.095M). As 
indicated, values computed at lower resolutions are rescaled such that they would coincide for 
perfect second order convergence. 



Note that a function g has been subtracted from the Brill wave function q. The form of go is 
found by setting the 3-metric equal to the spatial part of the Kerr metric and setting q — 0. 
Apart from the function go this is the same three-metric used for the non-rotating, distorted 
black hole spacetime. Here the logarithmic coordinate 77 is related to the Boyer-Lindquist 
coordinate r by 

r = r + cosh 2 (77/2) — r_ sinh 2 (??/2), (17) 

where r± = M ± V M 2 — a 2 . When considered together with the extrinsic curvature, this 
3-metric allows for a distinct class of initial data: distorted Kerr or Bowen and York black 
holes. To obtain initial data for a distorted Kerr black hole, one can set the conformal 
extrinsic curvature tensor to that which would be obtained for the standard Kerr metric, 
as shown in [Q. 

In order to obtain a distorted Bowen and York black hole, one can use the Bowen and 



York expression for the conformal extrinsic curvature [29 



3 J sin 2 6 
Kij=\ I , (18) 
3 J sin 2 9 
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where J is the total angular momentum of spacetimeS. The conformal metric used in this 
case is that resulting from setting q = in (16) which is related to the conformally flat 
metric. 

It is not possible to extend the axisymmetric distorted Bowen and York black holes to 
3D in the same way as in the non-rotating case. This is because if one allows the Brill wave 
function q to depend on the azimuthal angle <ft, it can be shown that the ^-component of 
the momentum constraint equation reduces to Jd^q = 0; i.e. for the Bowen and York form 
of the conformal extrinsic curvature, one cannot have both angular momentum and a Brill 
wave function with azimuthal dependence. Instead, to extend the rotating case to 3D, we 
use the conformally flat metric that results from setting q = q in equation (0), and place 
the 0-dependence in the conformal extrinsic curvature. The general form of @ used is 



Kij = m sin(m<; 



f h h 
hi h x -h 
\ -h x sm 2 9 



( o o k V(j> \ 
o o k H 

K v< j, Ke^ J 



(19) 



where 



K v 



sin 2 9[3J + cos(m( l 



1 



sin 3 ^ 

cos(m0) sm9(m 2 v — sin 2 9d v g). 



<%(sin 4 eg)}, 



(20) 
(21) 



Note that we have made all ^-dependence in the problem explicit. The functions v, g, ho, h x , 
and hi that appear above depend only t] and 9. Below we outline a procedure for finding 
forms for these functions that result in a conformal extrinsic curvature that satisfies the 
momentum constraint. The quantity J is a constant. 
We choose the functions h and hi to have the form 



h = d„Vt ^— -c^(sin0A), 



hi = d v A + 



sm9 



d e (smen). 



(22) 
(23) 



1 J is both a free parameter and the total angular momentum (ADM angular momentum) of the 
initial data, and hence of the space time. This is always true. It does not matter whether the 
data is axisymmetric or not, because the total angular momentum is defined at infinity (i.e. it is 
not local). Notice that when there are symmetries, it is possible to define not only the angular 
momentum at infinity, but also a quasi-local angular momentum by integrating the Killing vector 
(and both coincide). However, even this quasi-local definition is possible for our case: for the 
non-axisymmetric extrinsic curvature described in this section, since our data is conformally flat 
and a maximal slice, one can use the conformal Killing vector to define this quasi-local angular 
momentum as described in Ref. 

2 There is a way to construct "all" the solutions in an explicit form: namely construct all the 
solutions in terms of derivatives of two free functions as in Theorem 14 |3l|] 
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If m = 0, the momentum constraint equations are satisfied trivially. If m ^ 0, the 77, 0, and 
component of momentum constraints become, respectively, 

(Sg — 1 + $ + 2cot00fl)fi = (4cos0 + sin0<%)#, (24) 
(Sj - 1 + dg + 2 cot 0d fl ) A = - esc 8m 2 v - d e h x - 2 cot 8h x + sin 9d v g, (25) 
/i,, + 2 cos 8v + sin 0<%w = 0. (26) 

Now we expand all functions in trigonometric functions of 8 to the lowest order at which 
everything can be solved. After some calculations, we found the minimum collection of 
terms for g and Q to be 

9 = (90 + £ 2 sin 2 0)cos0 (27) 
Q = sin 2 8 tt 2 + sin 4 8 Q 4 (28) 

Note that for this discussion, tildes are used to indicate functions which depend on 77 alone. 

One substitutes these equations into the 77— component of the momentum constraint (PI). 
By doing this, one obtains an equation which can be written 

= E + E 2 cos (20) + £ 4 cos (40) . (29) 

Obviously, each coefficient function Ei must vanish. Setting E4 to zero gives 

g 2 = 1(25^- d 2 M. (30) 



Setting Eq + E 2 = gives 
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~g = l& 2 . (31) 



Next, we introduce an arbitrary function y, and redefine Q 2 to be 



fi 2 = y - j&4- (32) 



Using the relation E 2 — Eq = we obtain 

7 



Jl 4 = -(2^-3j?). (33) 



At this point, we have used the 77— component of the momentum constraint ( p4[ ) to define 
the functions Q and g in terms of the arbitrary function 2/(77). 

In order to solve the 0— component of the momentum constraint (p5|), we expand the 
functions A and v as follows: 

A = sin (20) A 2 + sin (40) A 4 , (34) 
v = {;[cos0-cos(30)] . (35) 

Plugging these expansions into the 0— component of the momentum constraint, and using 
the 0— component of the momentum constraint to write h x in terms of v, results in an 
equation which can be written 
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= L\ cos 9 + L 3 cos (39) + L 5 cos (50) . (36) 

Clearly, it must be the case that Li = L 3 = L 5 = 0. Setting L\ + L 3 + L 5 = gives 

A 2 = -2A 4 . (37) 

In order to solve for A 4 , we set L 1 /32 + L 3 /48 = 0, which eliminates the derivatives of A 4 . 
This gives 

~ 4(m 4 - I6)v + 15d v y - 9d?y 

A 4 • (38) 

The next step is to define v and y in terms of an arbitrary function x: 

v = a nd^x, (39) 

n=l,3,5 

y = x + a 7 d^x. (40) 

Finally, setting L 5 = and using the above substitutions produces an equation which can 
be written 

A 1 d v x + A 3 d 3 v x + A 5 d 5 rj x + A 7 d 7 v x = (41) 
Solving the linear system that results from setting each A{ to zero gives 

ai = — — — , (42) 

32 - 20m 2 ' V ' 

° 3 = 80 - 50m 2 ' (43) 
3 

° 5 = 20(-8 + 5m 2 )' (44) 

a 7 = 16 (45) 
7 -40 + 25m 2 v ; 

We now have analytic forms for the equations v, g, ho, hi, and h x which are functions of 
an arbitrary function x(r]) and which satisfy the momentum constraint equations. In this 
work, we choose the function x to be the Brill wave function without the sin n 9 factor: 

x = Q [e~ ( - > + e~ { - > }. (46) 

Note that if we set the parameter m to zero, the problem reduces to an axisymmetric 
problem, although not to the one studied by Brandt and Seidel. If we further set Qo = 0, 
we recover the Bowen and York solution. 

Using the conformal extrinsic curvature derived above, we can solve the Hamiltonian 
constraint numerically for the conformal factor ip: 

drf + 39* +C ° tU d9 +CSC % 2 " C 8 13 * 10 

We then interpolate these data onto a 3D Cartesian grid. 

Fig. shows the logarithm of the residual of the Hamiltonian constraint, H, on the 
a;— axis at three different resolutions for the data set (Q ,r] ,a, J,m) = (1,0,1,35,2). The 
values of H are rescaled so that they would coincide for perfect second order convergence. 
As in the non-rotating case, we see near second order convergence, especially for the higher 
two resolutions. 
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FIG. 2. We show the convergence of the residual of the Hamiltonian constraint for distorted, 
rotating black hole initial data with the parameter set Q$ = 1.0, r/o = 0,cr = 1, J = 35.0, and m = 2. 
# is plotted for Cartesian grid sizes (resolutions) 259 3 (0.024M), 131 3 (0.047M), and 67 3 (0.095M). 
As indicated, values computed at lower resolutions are rescaled so that they would coincide for 
perfect second order convergence. 
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III. A SURVEY OF DISTORTED BLACK HOLE INITIAL DATA SETS 



In this Section we do a parameter study of the above initial data. For all cases considered, 
we choose the wave location and width parameters 7] and a to be and 1, respectively. 
This corresponds to a Brill wave located on the throat with unit width. This choice is 
made because we are concerned about the necessary closeness of the outer computational 
boundary for the future nonlinear evolutions of these data sets. 



A. Analysis of the ADM mass 



Although one cannot define a local energy density of the gravitational field in general 
relativity, if the spacetime is asymptotically flat, one can define the energy of an isolated 
source as measured by a distant observer. Arnowitt, Deser, and Misner defined an energy in 
natural way for their 3 + 1 formalism | 3~2[| . O Murchandha and York modified their expression 
to use the variables in York's conformal decomposition method |3*3]] . For conformal metrics 
which fall off fast enough with radius, their expression for the ADM mass is 



M = 

2"7T joo 



m Joo 



(48) 



Since our numerical domains have finite extent, we perform this integral over a sphere of 
large constant radius R, resulting in the following expression: 



M = — 



R 2 r 27r r* dip 
2ix Jo Jo dr 

In the ^-coordinate system, this equation is 



sin0d0. 



(49) 



M = -- 



2vr 



I Mo 

„ Vmax 
2 JO 



71 fdifj y 
o Vdrj ~ 2 



sin6d6. 



(50) 



When we state the ADM mass for a spacetime, we have confirmed that the error involved in 
computing the integral at a finite radius is not significant by making sure that the integral 
as a function of R is approaching a constant value where we are performing the integration. 

For distorted non-rotating cases, we studied a variety of Brill wave amplitudes (Q ) and 
azimuthal factors (c). To better understand the results, we first note that the azimuthal 
factor in the Brill wave function can be written 



1 + c cos 



| + |cos(20). 



(51) 



That is, it has a constant part and an oscillatory part. The constant part by itself in effect 
increases the Brill wave amplitude. 

Fig. |3] shows plots of ADM masses as functions of Qo and c. In panel 1 (for c = —2), 
the constant term in the azimuthal factor vanishes. One expects that for this case, one will 
obtain the same spacetime if one changes the sign of the Brill wave amplitude, only rotated. 
Thus, the ADM mass should be symmetric about Q = 0, which is seen. In panel 2 (for 
c = —1), we see the results consistent with an effective amplitude decrease. For Q < 0, 
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0.6 -0.3 0.3 0.6 

(4)c=l 



-0.6 -0.3 0.3 0.6 -0.6 -0.3 0.3 0.6 

Qo 

FIG. 3. We show the ADM mass as a function of the Brill wave parameters Qo for four different 
values of c. The wave location and width parameters rjo and a are and 1, respectively. Grid size 
(resolution) is 131 3 (Ax = 0.2). 



the mass increases with decreasing Q , and has lower values than the corresponding Q in 
axisymmetry. For Qq > 0, the minimum in the mass occurs at a higher value of Qo than 
the axisymmetric case, and it has a lower value at any given c. In panel 3 (c = 0), the 
results are what one might expect given Bernstein's results in axisymmetry, and the form 
of the azimuthal factor in the g-function. Bernstein found that in axisymmetry, when one 
increases Qo, the ADM mass initially decreases, then increases with Q . Although the initial 
increase seems counter-intuitive, he points out that when one increases the amplitude of the 
Brill wave, it is possible to simultaneously decrease the mass of the underlying black hole, 
which can give a net decrease in the mass. Bernstein also found that for the negative Brill 
wave amplitudes, the ADM mass increases monotonically with decreasing Q . Our results 
for c = agree with the above Bernstein's results. Panel 4 (c = 1) shows the same features 
as the panel 3. However, the minimum in the ADM occurs at a lower value of Qo and is 
larger at Qo = 0.6. From Eq. fll4|), the case of Qo = corresponds to a Schwarschild black 
hole. For this case, the ADM mass should be the Schwarzschild mass Mo, which in this work 
is chosen to be 2. This is the value seen. 

The effect of the oscillatory part of the azimuthal factor by itself can be seen in panel 1 
of Fig. |3|. Increasing the absolute value of c increases the ADM mass. For positive values c, 
both the constant term and the oscillatory term are increased, and the effect of the ADM 
mass has the same character as increasing the Brill wave amplitude. For negative values of c, 
the situation is more complicated. For the range —2 < c < 0, the magnitude of the constant 
term decreases with decreasing c while the magnitude of the oscillatory term increase. The 
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FIG. 4. We show the ADM mass M as a function of the total angular momentum J for various 
various values of the Brill wave amplitude Qq. The solid lines show the axisymmetric cases (m = 0), 
and the dashed line shows the non-axisymmetric cases (m = 2). 

relative effect on the ADM mass within this range depends on the sign of the Brill wave 
amplitude. 

For distorted rotating cases, we investigate how the ADM mass behaves as a function 
of the Brill wave amplitude (Qq), the total angular momentum (J), and the azimuthal 
dependence variable (m). 

Fig. | shows the ADM mass M as a function of the total angular momentum J for 
various values of the Brill wave amplitude Qq. The ADM masses increase monotonically 
with increasing amplitude for both axisymmetric and non-axisymmetric cases. When the 
total angular momentum J is small, it is not effective for increasing the ADM mass. 

We also observe the different results when comparing axisymmetric and non- 
axisymmetric cases. In the axisymmetric cases (m = 0), the non- vanishing components 
of the conformal extrinsic curvature are K v< p ( p0|) and Kg^ ([21]) . By substituting Eq. ( |46| ) 
into Eqs. ([20D, we see that the component can be written 




0) = f + h cos (26) + U cos (40) + /„ cos (69) 



(52) 



where 




(55) 



(53) 



(54) 
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h = ^92 (56) 
On the other hand, for m = 2, this component takes the form 

K v4> {m = 2) = f Q + f[ cos (20) + f 2 cos(20 - 69) + f z cos(20 - 40) + f A cos(20 - 29) + 

# cos(20 + 26) - ft cos(20 + 49) + f 2 cos(20 + 69) (57) 

where 

# = \j (58) 

^ = S^-5 J (59) 

£ = lih (60) 

^ = Y6~ 90 - k 92 (61) 

/I = \h - 1,92 (62) 



If we compare the constant components of Eq. (p2|) and Eq. (0), we see that fo incorporates 
both the total angular momentum J and the amplitude Qq, whereas /q incorporates only J. 
Because of this, for low spin, the ADM mass increases more rapidly with amplitude for the 
axisymmetric cases, as can be seen in Fig. [|. As J increases, we see the effect of the Brill 
wave amplitude decreasing and approaching the pure Bowen-York values. 

Fig. | shows J/M 2 as a function of the total angular momentum J for various values of 
the Brill wave amplitude Q . Whereas the ADM mass increases more rapidly with amplitude 
for the axisymmetric cases, J/M 2 increases more rapidly for non-axisymmetric cases. Also, 
as J increases, we see the effect of the Brill wave amplitude decreasing and approaching 
the pure Bowen-York values. This is consistent with the above ADM mass results. Notice 
that for the Bowen-York data, J/M 2 does not converge to unity as discussed in This 
indicates that even for non-axisymetric cases, J/M 2 for our initial data will never reach unity. 
However, it is interesting to point out that if we add a non-axisymmetric perturbation to 
the second fundamental form in [34], it would cause J/M 2 to increase even more. 



B. Analysis of apparent horizons 

In this Section we discuss the location and properties of the apparent horizon (AH) in 
these initial data sets. Studying the apparent horizon gives us a detailed understanding of 
the shape and mass of the black holes. 

Defining s M to be the outward-pointing space-like unit normal of a two-sphere S embed- 
ded in a constant time slice S with timelike unit normal n M , we can construct the outgoing 
null normal to any point on S as k^ 1 = n M + s M . The surface S is called a marginally trapped 
surface if the divergence of the outgoing null vectors V»A; M vanishes, or equivalently, if 



= As' + K ij s i s j - K = 0, (63) 
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FIG. 5. We show J/M 2 as a function of the total angular momentum J for various various 
values of the Brill wave amplitude Qq. The solid lines show the axisymmetric cases (m = 0), and 
the dashed line shows the non- axisymmetric cases (m = 2). 



where G is the expansion of outgoing rays. The AH is defined as the outer-most marginally 
trapped surface. 

By construction, all of our data sets contain an apparent horizon. Gibbons showed that 
an isometry surface, such as the throat, is an extremal area surface [06 ]. This means that 



in the some cases (for example, the time symmetric cases) the throat will be an apparent 
horizon. But, in general, the apparent horizon will not be on the throat. 

We have computed the coordinate location of the apparent horizon in our initial data 



sets using the flow algorithm described in [37 



Fig. ^| shows the coordinate location of the apparent horizon for several non-rotating 
initial data sets. The parameters are (Q , c) = (Q , 1) with the following values of Q 
from the outer to the inner: Q = (—0.6, —0.5, —0.25,0.25,0.5). We observe that data sets 
with Qo = (—0.25, 0.25, 0.5) have AH's which are on the throat, while data sets with Q = 
(—0.6, —0.5) form new horizons. Although the above results are only for a few examples, 
we observe that the azimuthal factor is large enough that the effective amplitude increase 
results in another minimal surface forming. 

Fig. |7| shows the coordinate location of the apparent horizon for the axisymmetric rotating 
cases with parameters (Qo, J) = (1,J)- The solid (outermost) line corresponds to J = 0. 
When J is increased, the AH tends to shrink. For high spin cases, the horizon will shrink. 
This behavior is similar to the location of the AH of Kerr. Thus, the AH's should be on the 
throat for the above parameter sets. 

Fig. ^| shows the coordinate location of the apparent horizon for non-axisymmetric rotat- 
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FIG. 6. The coordinate location of the apparent horizon is plotted for amplitudes 
Qo = (—0.6,-0.5,-0.25,0.25,0.5). The other parameters are (rjo,cr, c) = (0,1,1). The outermost 
AH corresponds to Qo — —0.6, and the dashed line corresponds to Qo = —0.5. The remaining data 
sets have an AH at the inner circle, which is the throat. Grid size (resolution) is 131 3 (0.2). 



ing cases with parameters the same as in the axisymmetric case: (Qo, J) = (1, J). The dot 
and dashed line (outermost) corresponds to J = 75 (J/M 2 = 0.9). For the rapidly spinning 
cases, horizons are expected to expand. This means that compared with axisymmetric cases, 
when we increase the Brill wave amplitude, the AH horizons will leave the throat easily. 

When the apparent horizon is on the throat, it is a coordinate sphere. It is interesting, 
however, to study how distorted it is in the physical space. One measure of this distortion 
is the ratio of the polar to equator circumference: C r = Cp/Ce- The polar circumference, 
Cp, is the proper length of the circumference which goes through the poles along a line of 
constant 0, and the equatorial circumference, Ce, is the proper length of equator, denned 
as 9 = 7r/2. In axisymmetry, this ratio is independent of <ft, and is a single number which 
characterizes a data set once one has established the above definitions. In this case, if 
Cp/Ce = 1, the horizon is spherical, if Cp/Ce < 1 it is oblate, and if Cp/Ce > 1 it 
is prolate. In 3D, one can establish analogous definitions, although in this case Cp, and 
therefore C r , is a function of <fi. 

Fig. |9] shows the ratio Cp/Ce for two values of <p for non- rotating cases with parameters 
(Qo, Vo, 0"; c) = (Qo, 0, 1, c). Cpi is the value of Cp at = 0, and Cp 2 is the value of Cp at 
(ft = 7r/2. We observe that highly distorted features appear when c > and Qo > 0. 

For distorted rotating black holes, we study the same parameter sets as above: 



(Qo, ^0) 0, J, m ) = (1, 0, 1, J, m). Fig. [10] shows C r as a function of J for m = and m — 2. 



Although the distortion patterns are different for axisymmetric and non-axisymmetric cases, 
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FIG. 7. The coordinate location of the apparent horizon of axisymmetric cases for various total 
angular momentum parameter J :150, 30, and which correspond J/M 2 : 0.9, 0.5, and respectively. 
The solid line (outermost) is J = 0. When J is increased, the horizon shrinks. 

both cases are nearly spherical for the range of parameters studied. High spin axisymmetric 
cases are slightly more oblate than the other cases. Also, for the non-axisymmetric cases, 
Cpi and Cp2 are indistinguishable. 



C. Maximum Radiation Loss (MRL) 



Another quantity of interest which is related to the apparent horizon is the maximum 
radiation loss (MRL). The second law of black hole thermodynamics states that the area of 
the event horizon of a black hole, Aeh-, cannot decrease in time. One can define a quantity 
associated with the area of black hole known as the irreducible mass, M irr 



Mi, 



I A 



EH 



16tt 



(64) 



By analogy to this quantity, one can define the mass of the apparent horizon given its area, 

A AH, 



M 



'A 



AH 



AH 



167T 



(65) 



More generally, for a rotating black hole, one can define the apparent horizon mass as 
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FIG. 8. The coordinate location of the apparent horizon of non-axisymmetric cases for various 
total momentum parameter J: 75, 10, and which correspond J/M 2 : 0.9,0.5, and respectively. 
Contradictory to axisymmetry cases, the outermost one is high spin: J = 75. 



, A 2 A AH . 4ttJ 2 

M AH = -T7— + ~ A , (66) 

167T A AH 

where J is the total angular momentum. 

One defines the radiation efficiency RE as the difference between the ADM mass of a 
spacetime and its irreducible mass in its final state: 

M - M irr . f , . 

RE = (67) 

However, by second law of black hole thermodynamics, the final irreducible mass must be 
greater than or equal to the initial irreducible mass. Thus, 

RE < M ~ A f-. (68) 

Because one needs to perform an evolution to find the location of the event horizon on any 
slice, we can not compute the above fraction, even though it is for quantities defined on 
our data set. However, for our initial data sets, we know that the apparent horizon is the 
outer-most minimal area surface. It can be shown that if an apparent horizon exists, it must 
lie inside or coincident with an event horizon j|C]. Thus the area of the event horizon on 



our initial slice must be greater than or equal to the area of the apparent horizon. Then we 
can define a quantity we call the maximum radiation loss (MRL) by using the ADM mass 
M: 
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FIG. 9. The ratio Cp/Ce of non-rotating case for the parameter sets (Qo, Vo, cr, c) = (Qo, 0, 1, c). 
Cpi is the value of Cp at <fi = 0, and Cp2 is the value of Cp at cj> = tt/2. 



MRL 



M - M AH , 
M 



> RE. 



(69) 



The fraction above is an upper bound on the radiation efficiency. Therefore, we know that 
during an actual evolution, the total energy radiated has to be less than this amount. This 
will be a guide to choosing interesting initial data sets for evolution in the future. 



Fig. 11 shows the MRL for non- rotating function of Qo and c for the family 

of initial data sets (Qo,r]o,a,c) = (Qo,0, l,c). The largest MRL is about 10% of the total 
mass. The sets within this family with the largest MRLs are those with both positive Qq 
and c, indicating that these might be interesting data sets to evolve. 

Fig. O shows the MRL for rotating CcLSGS cLS cL function of J and Qo f° r the family of 
initial data sets (Q , rj , a, J) = (Q ,0,1, J). For the axisymmetric cases, when the initial 
amplitude is larger the MRL tends to decrease, while for non-axisymmetric cases the trend 
is reversed for most values of J. Also, for the axisymmetric cases, the MRL increases as 
J becomes larger, but for non-axisymmetric cases the MRL decreases. As we saw for the 
ADM masses, when J is increased, the MRL for both classes of initial data approaches the 
pure Bowen and York values. 

Furthermore, we can address an important question. Penrose proposed a criterion for a 



spacetime which, if violated, would indicate that Cosmic Censorship would be violated [41 



The criterion, the so-called Penrose inequality, is simply that the irreducible apparent hori- 
zon mass must be less than or equal to the ADM mass. This would result in a negative value 
of the MRL. Using the Brill wave plus black hole initial data, violations of the Penrose 
inequality were looked for in Ref . [P]fl[|2J , but none were found. No violations were observed 
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FIG. 10. The ratio Cp/Ce for axi- and non-axisymmetric rotating cases. The parameter sets 
are (Qo, rjo, o, J) = (1, 0, 1, J). For non-axisymmetric cases, we have Cp\ which is the value of Cp 
at <j) = 0, and Cpi which is the value of Cp at (j) = tt/2. The solid line shows the axisymmetric 
case, and dashed lines show non-axisymmetric cases. High spin axisymmetric cases are slightly 
more oblate than the other cases. 

in either our non-rotating or rotating initial data sets. 



IV. CONCLUSIONS 

We have extended Bernsteins's axisymmetric, distorted black hole data sets to full 3D 
by giving them a particular azimuthal dependence in the Brill wave function q. For rotating 
initial data, we have generalized axisymmetric distorted rotating black hole by allowing the 
conformal extrinsic curvature to have an azimuthal dependence. 

From their physical properties, such as their ADM mass distortion of horizons and MRL, 
non-rotating axisymmetric cases are consistent with what has been observed by Bernstein's. 
In particular, Bernstein's results are recovered when the azimuthal parameter c is set to zero. 
However, rotating data sets are, by construction, different from those of Brandt and Seidel. 
For those initial data sets, we have seen that one of distortion parameters, Qo, strongly 



effects the physical properties when J is low, as we showed in Section III. However, as J 



is increased, physical properties are approaching the pure Bowen-York data. For further 



investigation of higher amplitude and spin cases for full 3D initial data as in Ref. pi 
since we are using uni-grid Cartesian coordinate for all physical quantities except for the 
ADM mass, we need larger grid size and higher resolutions by using such as adaptive mesh 
refinement. 
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FIG. 11. We show the MRL for the initial data set: {Qo,1]o,<t,c) = (Qq,0,1,c). The efficient 
MRLs are for the positive amplitude and azimuthal factor. 



In future papers we would like to address the radiation of angular momentum by evolving 
non-axisymmetric, distorted, rapidly-rotating black holes. Direct comparisons of radiative 
energy were done for axisymmtric and non-rotating initial data [p0 , pl[ 22,43 . The close limit 
approximation studies are addressing radiated angular momentum [ 44|, 45 1 . However, as men- 
tioned in Ref. a direct comparison with full numerical simulations is under investigation. 
By using the Lazarus method [||(|,(|7j, Ref. |L5|, one shows that initial data for grazing col- 
lisions cannot be mapped into the perturbative method. Although our non-axisymmetric, 
distorted rapidly rotating black holes are not necessarily astrophysically relevant, such anal- 
ysis will provide an example of the usefulness of perturbation theory as an interpretive tool 
for understanding the dynamics produced in full nonlinear evolutions. 
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FIG. 12. We show the MRL for rotating cases. The initial data sets are 
(Qo, tjq, a, J) = (Qo,0, 1, J). The solid lines show the axisymmetry cases, and the dashed line 
shows the non-axisymmetry cases. 
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